Phot on- number correlations by photon-number 

resolving detectors 



Alessia Allevi, 1 Maria Bondani, 2 '* and Alessandra Andreoni 3 

1 C.N.I.S.M., U.d.R. Milano, 1-20133, Italy 

2 National Laboratory for Ultrafast and Ultraintense Optical Science - C.N.R.-I.N.F.M., 

Como, 1-22100, Italy 

3 Dipartimento di Fisica e Matematica, Universita degli Studi dell'Insubria and 
C.N.I.S.M., U.d.R. Como, 1-22100, Italy 

* Corresponding author: maria.bondani@uninsubria.it 



We demonstrate that by using a pair of photodetectors endowed with 
internal gain we are able to quantify the correlation coefficient between the 
two components of a pulsed bipartite state in the mesoscopic intensity regime 
(less than 100 mean photons). © 2010 Optical Society of America 
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The existence of photon-number correlations is necessary to produce conditional states 
based on number photodetection. Moreover the study of the properties of photon-number 
correlations and, in particular, of the distribution of the difference in the number of photons 
detected in the two components of a pulsed bipartite state, is meaningful of the properties 
of the state [1,2]. For instance, measurements of sub-shot noise correlations have been used 
to determine the nonclassicality of both bipartite [3-5] and tripartite states [6] . 
Different types of commercial detectors can be used to measure photon numbers in the differ- 
ent intensity regimes: pin photodiodes [1,7], amplified pin photodiodes [3,5,6], CCD's [4,8,9], 
loop detectors [10,11], hybrid photodetectors [12-14], multipixel photodetectors [15,16]. Some 
of them have been exploited to reconstruct photon-number statistics [3,6,8,11-16] or for cor- 
relation measurements [3-6, 8,9]. 

We have already demonstrated that hybrid photodiode modules (HPD, R10467U-40, Hama- 
matsu, maximum quantum efficiency ~ 0.4 at 550 nm) and Silicon photomultipliers (SiPM, 
MPPC S10362-11-100C, Hamamatsu, maximum quantum efficiency ~ 0.65 at 440 nm) are 



endowed with the capability of reconstructing the correct statistics of detected photons pro- 
vided their outputs are analyzed in the proper way [12-14,16]. In this Letter we exploit the 
fact that, in principle, our analysis allows us to determine the number of photons detected 
at each laser shot and demonstrate that pairs of the above detectors can be used to measure 
shot-by-shot photon-number correlations between the components of a bipartite state. To 
test the detection procedure, we start by measuring the correlations in a classically-correlated 
bipartite state, namely a single-mode thermal field divided by a beam-splitter. 
The experimental setup is sketched in Fig. 1. The single- mode thermal field is obtained by 
passing the second-harmonic pulses (~5.4 ps, 523 nm) of a Nd:YLF mode-locked laser am- 
plified at 500 Hz (High Q Laser Production) through a rotating ground-glass plate (P in 
Fig. 1) [17] and selecting a single coherence area in the far field by a pin-hole (PH, 150 /im 
diameter). The selected mode is then sent to a polarizing beam-splitter (PBS) whose outputs 
are collected by two multimode fibers and delivered to the detectors. To obtain a fine tuning 
of the beam-splitter transmittance we inserted a half-wave plate (HWP) between the PH 
and the PBS. Single-shot detector outputs were then amplified, integrated, digitized and 
recorded. We measured 30000 subsequent shots with the HPD detectors and 50000 with the 
SiPM detectors at a number of different intensity values set by a variable filter (F). Figure 2 
shows two typical pulse height spectra obtained at similar intensities: clearly, the response 
of the SiPMs produces better resolved peaks as compared to the HPDs. As the HPDs are 
not affected by dark-counts, by applying the selfconsistent procedure described in [12] to the 
output of the detectors, we can evaluate the number of detected photons in each laser shot. 
The procedure gives the value of the overall gain (7) of the detection apparatus that links 
the final outputs (voltages for HPDs and multichannel addresses for SiPM) to the number 
of detected photons. In the case of SiPMs, the procedure described in [16] yields a value of 
detected photons that includes the unavoidable effects of dark-counts and cross-talk, that 
must be taken into account to properly interpret the experimental results. The procedure 
itself yields the values of the dark-count rate Wide and of the cross-talk probability e for the 
considered data-set. 

In Fig. 3 we plot the joint probability distributions obtained by re-scaling the outputs of 
the detectors by the proper value of 7 but without rebinning the values. As expected, the 
results for SiPM display a better localization of the peaks, which should promise a better 
determination of the correlation. 

Referring to Fig. 1, the field operators at the outputs of the beam splitter (c and d) are 
linked to the input ones (a and b) by: 

c = y/ra — a/1 — t b , d = y/1 — r a + y/rb . (1) 

If b is in the vacuum, the moments of the bipartite state at the output are given by (n c ) = 
r(n a ), < = r [r< + (1 - r)(n a )], (n d ) = (1 - r)(n ), < = (1 - r) [(1 - r)< + T (n a )} 
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and (n c n d ) = t(1 — r) ((n 2 ) — (^a))> from which we obtain the expression for the normalized 
correlation coefficient: 

th _ ((n c ~ (n c )) (n d - (h d ))) 



<T 2 <7 2 



y/frvl + (1 - r)(n a >] [(1 - r) < + r<n a >] 



(2) 



For a single-mode thermal state at the input, p a = J2 n P^( n )\ n ) n ( n \ with p t h(n) = W 1 /(l + 
rl) 1+n and (n a ) = n, we have 



V(l - r)n 2 \Jn c n d 



v/(rn + l)[(l-r)n + f] ^ + f )(7T d + 1) ' 

where n c = (n c ) and n d = (n d ) . If the primary detection event is performed with quantum 
efficiency rj and no dark counts, the expression in Eq. (3) remains the same upon substituting 
the mean number of photons (n c ,d) with that of detected photons (m c d = r] c ^ d n Cjd ) [1] . 
Figure 4 shows the results for shot-by-shot correlation measured with HPDs and SiPM (dots) 
as a function of the mean value of output detected light (m c + m d ) . The full lines are the 
theoretical values evaluated according to Eq. (3) by inserting the experimental mean values. 
Not surprisingly the results obtained with HPDs perfectly superimpose to the theory, while 
those obtained with SiPMs display a reduced correlation, due to the effect of dark-counts 
and cross talk. By evaluating the theoretical curve at the effective mean values m e ff = (m — 
™dc) /(I + e), we obtain the dashed curve that fits the experimental data rather satisfactorily. 
Relevant applications of correlated bipartite states require the determination of the nature 
of the correlations. A quite standard way to discriminate between classical and nonclassical 
states [1] is to study the statistics of the difference, p(S), in the number of photons detected 
shot-by-shot at the outputs of the beam splitter. Figure 5 displays the experimental results 
for p(S) evaluated by subtracting the number of detected photons measured shot-by-shot 
with HPDs and SiPMs together with the theoretical values evaluated according to Ref. [1]. 
The sample curves where chosen to have very similar detected mean values for the two classes 
of detectors. By evaluating the fidelity (/ = J2i \/Pi(^)Pi h (^)) we n °t e that quality of the 
results is almost the same for HPDs and SiPMs, although for the first is slightly better (see 
figure caption for the values of /). The Insets of Figure 5 display the measured variance 
(T 2 (5) as a function of the total number of detected photons, as compared to the theoretical 
predictions u 2 {5) = (m c — m d ) 2 + m c + m d [1]. The correct evaluation of & 2 (5) is crucial for the 
estimation of the noise reduction in a bipartite state R = a 2 (S)/(fn c + fn d ), which is a means 
to discriminate classical from nonclassical states. Again the results obtained with HPDs 
superimpose to theory, whereas those obtained with SiPM follow a broader distribution: 
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these results indicate that HPDs can give direct information on sub-shot noise level, whereas 
the information obtained with SiPM must be corrected to be interpreted. 
In conclusion we have demonstrated that detectors partially endowed with photon-number 
resolution can be used to measure photon-number correlations shot-by-shot and to reliably 
estimate the noise reduction factor. We have also demonstrated that, although endowed with 
a better photon-number resolution, SiPM are affected by spurious effects (dark-counts and 
cross-talk) that degrade the quality of results and require corrections that risk to prevent 
their use for state preparation by conditional measurements [15]. However, with both kinds 
of detectors our system is less cumbersome than that described in Ref. [10] and, in the case 
of HPD, provides similar overall detection efficiency 
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discussions. 

Present address of Alessia Allevi is C.N.I.S.M., U.d.R. Como, 1-22100, Italy. 
References 

1. A, Agliati, M. Bondani, A. Andreoni, G. De Cillis, and M. G. A. Paris, "Quantum and 
classical correlations of intense beams of light investigated via joint photodetection" , J. 
Opt. B: Quantum Semiclass. Opt. 7, S652S663 (2005). 

2. I. P. Degiovanni, M. Genovese, V. Schettini, M. Bondani, A. Andreoni, and M.G.A. Paris, 
"Monitoring the quantum-classical transition in thermally seeded parametric down- 
conversion by intensity measurements", Phys. Rev. A 79, 063839 (2009). 

3. M. Bondani, A. Allevi, G. Zambra, M. G. A. Paris, and A. Andreoni, "Sub-shot-noise 
photon-number correlation in a mesoscopic twin beam of light" , Phys. Rev. A 76, 013833 
(2007). 

4. O. Jedrkiewicz, Y. K. Jiang, E. Brambilla, A. Gatti, M. Bache, L. A. Lugiato, and P. Di 
Trapani, "Detection of sub-shot-noise spatial correlation in high-gain parametric down 
conversion", Phys. Rev. Lett. 93, 243601 (2004). 

5. I. N. Agafonov, M. V. Chekhova, and G. Leuchs, "Two-Color Bright Squeezed Vacuum", 
quant-ph/0910.4831v2 

6. A. Allevi, M. Bondani, M.G.A. Paris, and A. Andreoni, "Demonstration of a bright and 
compact source of tripartite nonclassical light", Phys. Rev. A 78, 063801 (2008). 

7. O. Aytiir and P. Kumar, "Pulsed twin beams of light", Phys. Rev. Lett. 65, 1551-1554 
(1990). 

8. O. Haderka, J. Perina, Jr., M. Hamar, and J. Perina, "Direct measurement and recon- 
struction of nonclassical features of twin beams generated in spontaneous parametric 



1 



down-conversion", Phys. Rev. A 71, 033815 (2005). 
9. G. Brida, L. Caspani, A. Gatti, M. Genovese, A. Meda, and I. Ruo Berchera, "Measure- 
ment of sub-shot-noise spatial correlations without background bubtraction" , Phys. Rev. 
Lett. 102, 213602 (2009). 

10. M. Micuda, O. Haderka, and M. Jezek, "High-efficiency photon-number-resolving mul- 
tichannel detector", Phys. Rev. A 78, 025804 (2008). 

11. W. Mauerer, M. Avenhaus, W. Helwig, and Ch. Silberhorn, "How Colors Influence Num- 
bers: Photon Statistics of Parametric Downconversion" , quant-ph/0812.3597vl 

12. M. Bondani, A. Allevi, A. Agliati, and A. Andreoni, "Self-consistent characterization of 
light statistics," J. Mod. Opt. 56, 226-231 (2009). 

13. M. Bondani, A. Allevi, and A. Andreoni, "Light statistics by non-calibrated linear pho- 
todetectors," Adv. Sci. Lett. 2, 463-468 (2009). 

14. M. Bondani, A. Allevi, and A. Andreoni, "Wigner function of pulsed fields by direct 
detection", Opt. Letters 34, 1444-1446 (2009). 

15. I. Afek, A. Natan, O. Ambar, and Y. Silberberg, "Quantum state measurements using 
multipixel photon detectors," Phys. Rev. A 79, 043830 (2009). 

16. M. Ramilli, A. Allevi, V. Chmill, M. Bondani, M. Caccia, and A. Andreoni, "Photon- 
number statistics with Silicon photomultipliers" , quant-ph/0910.4786v2 

17. F. T. Arecchi, "Measurement of the Statistical Distribution of Gaussian and Laser 
Sources", Phys. Rev. Lett. 15, 912-916 (1965). 



5 



List of Figure Captions 

Fig. 1. (Color online) Scheme of the experimental setup. P, rotating diffuser plate; PH, pin- 
hole; F, variable neutral density filter; HWP, half wave plate; PBS, polarizing beam splitter; 
L, lenses; MF, multimode optical fibers; D c and D d , photo-detectors. 

Fig. 2. (Color online) Pulse height spectra at the outputs of the beam splitter as measured 
by HPD and SiPM photodetectors at very similar intensities. 

Fig. 3. (Color online) Experimental joint probability distribution. Fig. 4. (Color online) 
Shot-by-shot correlation, T, as a function of the mean value of output detected light (dots) 
together with the theoretical curves (full lines). The dashed line in the right panel is the 
theoretical prediction corrected for dark counts (m^c — 0.071, = 0.034) and cross talk 
(e c = 0.10, e d = 0.08). 

Fig. 5. (Color online) Experimental results for p(5) (full circles) and theoretical predictions 
(full lines). The fidelities of the data are / = 0.9999, 0.9997, 0.9996, 0.9994, 0.9991, for HPDs 
and / = 0.9997,0.9994,0.9994,0.9993,0.9991, for SiPM. Insets: measured variance a 2 (5) of 
p(S) as a function of the total number of detected photons (dots) and theoretical predictions 
(full line). 
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Fig. 1. (Color online) Scheme of the experimental setup. P, rotating diffuser plate; PH, pin- 
hole; F, variable neutral density filter; HWP, half wave plate; PBS, polarizing beam splitter; 
L, lenses; MF, multimode optical fibers; D c and D^, photo-detectors. 
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Fig. 2. (Color online) Pulse height spectra at the outputs of the beam splitter as measured 
by HPD and SiPM photodetectors at very similar intensities. 
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Fig. 3. (Color online) Experimental joint probability distribution. 
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Fig. 4. (Color online) Shot-by-shot correlation, T, as a function of the mean value of output 
detected light (dots) together with the theoretical curves (full lines). The dashed line in the 
right panel is the theoretical prediction corrected for dark counts {m dcfi = 0.071, m dCjd = 
0.034) and cross talk (e c = 0.10, e d = 0.08). 
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Fig. 5. (Color online) Experimental results for p(S) (full circles) and theoretical predictions 
(full lines). The fidelities of the data are / = 0.9999, 0.9997, 0.9996, 0.9994, 0.9991, for HPDs 
and / = 0.9997,0.9994,0.9994,0.9993,0.9991, for SiPM. Insets: measured variance a 2 {5) of 
p(S) as a function of the total number of detected photons (dots) and theoretical predictions 
(full line). 
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